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ABSTRACT 

Turbulence in the solar wind is believed to generate an energy cascade that is supported primarily by Alfven 
waves or Alfvenic fluctuations at MHD scales and by kinetic Alfven waves (KAWs) at kinetic scales k^pt > 
1. Linear Landau damping of KAWs increases with increasing wavenumber and at some point the damping 
becomes so strong that the energy cascade is completely dissipated. A model of the energy cascade process that 
includes the effects of linear collisionless damping of KAWs and the associated compounding of this damping 
throughout the cascade process is used to determine the wavenumber where the energy cascade terminates. It is 
found that this wavenumber occurs approximately when \ j/Ld\ ~ 0.25, where Ld{k) and 'f{k) are, respectively, 
the real frequency and damping rate of KAWs and the ratio j/oj is evaluated in the limit as the propagation 
angle approaches 90 degrees relative to the direction of the mean magnetic field. For plasma parameters typical 
of high-speed solar wind streams at 1 AU, the model suggests that the KAW cascade in the solar wind is 
almost completely dissipated before reaching the wavenumber k±pi ~ 25. Consequently, an energy cascade 
consisting solely of KAWs cannot reach scales on the order of the electron gyro-radius, k±pe ^ 1. This 
conclusion has important ramifications for the interpretation of solar wind magnetic field measurements. It 
implies that power-law spectra in the regime of electron scales must be supported by wave modes other than 
the KAW. 

Subject headings: Solar wind — turbulence, magnetohydrodynamics, kinetic theory 



1. INTRODUCTION 



Rece nt pa pers by [S chekochihin et al.l (120091) . iHowes et alJ 
(l2008ah . and lSchekochihin et alJ (120081) have described a sce- 
nario for turbulence in collisionless magnetized plasmas, such 
as the solar wind, that can be briefly described as follows. 
The turbulence at large scales, scales larger than the ion in- 
ertial length c/cupi and the thermal ion gyro-radius pi, con- 
sists of an energetically dominant Alfven wave cascade that 
transfers energy from large to small scales. At scales on the 
order of the proton gyro-radius, k±pi ~ 1, the wavevector 
spectrum of the turbulence is highly anisotropic with energy 
concentrated in wavevectors nearly perpendicular to the mean 
magnetic field Bq so that fcj^ » fc|| . At fcj^pi ^ 1, two things 
happen. On the one hand, there may be some nonlinear effects 
that are not well understood. On the other hand, a significant 
fraction of the energy in the Alfven wave cascade excites a 
kinetic Alfven wave cascade (KAW cascade) that carries the 
energy down to scales on the order of the thermal electron 
gyro-radius where the turbulence is finally dissipated by col- 
lisionl ess Landau damp ing. 

Sah raoui et alJ (120091) used the above scenario to interpret 
spacecraft measurements of solar wind turbulence in the ki- 
netic range of scales k±pi > 1. These authors used the lin- 
ear KAW dispersion relation and damping rates to argue that 
the KAW cascade should reach electron scales before being 
strongly damped, noting that the relative damping rate j/lu 
does not become of order unity until k±pe ^ 1, where p^ is 
the thermal electron gyroradius. The purpose of the present 
paper is to point out that this argument is incomplete because 
it does not take into account the compounding of the damp- 
ing throughout the course of the cascade process. KAWs 
are damped by collisionless Landau and transit-time damp- 
ing throughout the entire wavenumber range of their existence 
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from k±pi ~ 1 to k±pi ^ 1. If the energy cascade process 
is thought of as taking place in a sequence of discrete steps, 
then there is dissipation at each step in the sequence and the 
effects of damping are compounded with each step, analogous 
to the way compound interest works. When this compound- 
ing is taken into account it is found that for typical solar wind 
plasma parameters near 1 AU the KAW cascade is dissipated 
before reaching electron scales. This conclusion has impor- 
tant ramifications for the interpretation of solar wind power 
spectra in the kinetic regime. 

In astrophysical plasmas, the effects of collisionless 
damping on the turbulence spectrum have been stud- 
ied primarily using quasilinear theory (see, for example, 
iMelrosd 11994*) and phenomenological models in which 
the energy cascade process is modeled as a diffusion pro 
cess in wave n umber space (for example. iLi et al.l 
'Stawicki et al.' '200 1'; 'Cranmer & va n Ballegooijenl 
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Jiang et al. 2009; Matthaeus et al. 2009). Oth er app roaches 
include weak turbulence theory (lYoonl 12006 1, 2007; Gal tied 
2006; Chandran 200 8) and gyro -kinetic theory (Howes et al] 
2008b; iSchekochihin et alJl2009l) . Here we adopt a somewhat 
simpler approach that, like the diffusion models, is based 
on an equation expressing the conservation of energy in 
wavenumber space. Generally speaking, the objective of 
any of these models is to describe the essential physics as 
accurately and economically as possible. 

In this study, the effects of damping on the KAW cascade 
are modeled using two complimentary approaches. The first 
approach is heuristic and shows how this damping takes place 
through a sequence of steps in wavenumber space whereby 
energy is damped at each step before being transferred to 
higher wavenumbers. The second approach is based on an 
equation for the conservation of energy in wavenumber space 
and two different expressions for the energy cascade rate that 
are similar to Kolmogorov's relation e = kE(k)/T. This 
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approach h as much in common with the cascade model de- 
veloped by iHowes et alj ( l2008ah and allows both the energy 
cascade rate e and the energy spectrum E{k) to be computed 
as functions of wavenumber throughout the inertial range and 
dissipation range. In this study, the main objective is to de- 
termine the point in wavenumber space where the KAW cas- 
cade terminates and, therefore, the energy cascade rate e is the 
quantity of primary interest. Assuming that the propagation 
angles of the waves are nearly perpendicular to Bq, analytic 
solutions for the energy cascade rate e may be derived that 
are convenient for purposes of analysis and prediction. These 
analytic solutions are new and are presented here for the first 
time. 

One of the principal assumptions of this work is that the 
damping rates of linear wave theory are applicable at kinetic 
scales, even though nonlinear interactions are still strong. Al- 
though this same assumption has been made before by many 
investigators (Leamo net al. 1998; Ouataert 1998; Garv 1999; 
Leamon et all 119991; lO uataert & Gruzinov 1999; Li et al. 
IStaw icki et al.l 12 001.; .Cranmer & van BallegooiienI 
iHow es et al. 2008^, there is no well estabUshed 
criteria for its validity. This is, perhaps, the greatest source 
of uncertainty for the theory presented here. It is also easy 
to envision a scenario in which turbulence at MHD scales is 
dissipated via collisionless magnetic reconnection involving 
structures (reconnection sites) that span lengthsca les from the 
ion in ertial length to the electron inertial length ( Birn & Priest! 
120071 Chapter 3). To some extent, this is a possible alterna- 
tive to the wave ca scade and damping scenario described by 
iSchekochihin et all 12009 ) and others. 

For the moment, adopting the KAW cascade and linear 
damping scenario and assuming that the wave damping rates 
of the linear Vlasov-Maxwell theory are valid in the kinetic 
regime, then for typical solar wind conditions at 1 AU the 
models derived here predict that the KAW cascade terminates 
at wavenumbers less than or equal to k±pi ~ 25 which im- 
plies that the KAW cascade cannot reach electron scales in 
the solar wind at 1 AU. These results agree with and are sup- 
po rted by th ose obtained previously using the cascade model 
of iHowes et al. (2008a) which show a transition to an ex- 
ponentially decaying magnetic energy spectrum that decays 
rapidly at around the same wavenumber. On the contrary, 
th e magnetic e nergy spectra in the gyro-kinetic simulations 
of lHowes et al.l (l2008 b) do not appear to show any deviations 
from power- law behavior in the range 1 k±pi <S there 
should be if kinetic damping is having a significant effect. 
However, because these scales are near the smallest spatial 
scales in the simulation where an artificial electron hypercol- 
lisionality takes effect, the simulation may not accurately de- 
scribe the physics when fc^pi ^ 8. 

One note on terminology. The term "energy cascade" or 
"cascade," for short, refers to a nonlinear energy transfer pro- 
cess that transfers energy from large to small scales or vice 
versa. This term usually refers to the inertial range where the 
effects of dissipation are negligible. Here, this term is also 
used to describe the nonlinear energy transfer in the dissi- 
pation range. This slight abuse of notation should cause no 
confusion. 

2. SIMPLE HEURISTIC MODEL 

We begin with some general considerations that apply to 
the entire paper Consider a homogeneous proton-electron 
plasma with a constant uniform background magnetic field. 
The equilibrium state is charge neutral, free of macroscopic 
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Fig. 1 . — Damping rate of the wave amplitude in one wave period a{k) = 
exp[27r7(A;)/tj(A;)]. Tlie propagation angle 9 = 84° is shown in blue and 
9 = 89.9° is shown in red. 



electric fields, free of electric currents, and characterized 
by isotropic Maxwell distribution functions for both parti- 
cle species. It is assumed that the wave amplitudes at ki- 
netic scales are small enough that linear wave theory ade- 
quately describes the collisionless damping process. For a 
given wavenumber k (real-valued), the real and imaginary 
parts of the wave frequency oj + ij may be computed nu- 
merically using the hot plasma dispersion relation (Stix 1992] 
chap 10). For k± ^ k^^, the KAW is uniquely identified by 
means of its asymptotic dispersion relation lu — k\\VA in the 
limit as fc^ — > with kj_/k\\ or the angle of propagation held 
constant. 

In general, the damping is minimized for all wavenumbers 
in the range from k±pi ~ 1 to k±pi ~ 40 when the angle 
of propagation 6 is close to 90°, that is, near perpendicular to 
Bq. It is important to note that for all angles in a sufficiently 
small neighborhood of 90°, the ratio of the damping rate to 
the real frequency, 7/aj, is approximately independent of the 
angle 9 so that all propagation angles in this range yield the 
same damping per wave period. In practice, this means that 
the damping per wave period cannot be reduced further by 
going from 89.9° to 89.99°, for example. For simplicity, it 
will be assumed in this section that 6 is close, but not equal 
to 90°. This range of near-perpendicular angles is physically 
relevant because as a consequence of three-wave interactions 
the energy cascade process creates a spectrum in which the in- 
equality k± ^ fc|| is increasingly well satisfied as the cascade 
progresses to higher wavenumbers. 

Investigation of the damping of the KAW cascade begins 
with the plot in Figure- 1 which shows the attenuation of the 
wave amplitude in one wave period, that is, exp(27r7/cLi) ver- 
sus k± for a fixed angle of propagation close to 90°. As 
just mentioned, the curve for 89.9° in Figure-1 applies for all 
angles sufficiently close to 90°. Now recall the criti cal bal- 
ance hypothesis of iGoldreich & Sridhaii (119951 119971) which 
states that in the inertial range the cascade time is equal to 
the Alfven wave period. Using this as a guide, we assume 
that the cascade time in the kinetic regime is equal to the wave 
period of the KAWs. Consider the sequence of wavenumbers 
fc„ — 2''^ko, where kopi = 1 and n = 0, 1, 2, . . .. In one cas- 
cade time, one wave period, the energy at scale /c„ cascades 
to scale kn+i but in the process that energy is damped by a 
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TABLE 1 
Cascade rates with 
compounding taken into 
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0.95 
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1 


0.87 
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2 


0.84 
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0.73 


48% 
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0.52 


25% 
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16 


< 0.4 


7% 
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32 




< 1% 



factor af^, where a„ = exp[27r7(fc„)/w(fc„)] is the attenua- 
tion of the wave amphtude in one wave period. If e(/e„) is 
the energy cascade rate at wavenumber fc„, that is, the rate at 
which energy reaches scale fc„, then the rate at which energy 
reaches the next scale fc„+i is 



--n+l 



(1) 



where £„ — £(fc„). This is the change in the energy cascade 
rate caused by linear damping after just one step in the cas- 
cade process. Starting from n = 0, the energy cascade rate 
after n steps is equal to 
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(2) 



This describes how the energy of the waves is damped as it 
cascades through the kinetic range of scales. To compute the 
wavenumber dependence of e{k) for the KAW cascade it is 
only necessary to compute the relative damping rate from 
the hot plasma dispersion relation and then use equation (|2]l. 

By way of illustration, consider the damping rates for 89.9° 
shown in Figure- 1. For each wavenumber fc„ — 2"fco it is 
straightforward to read off the values of a„ from the plot in 
Figure- 1 and then compute the cascade rate e„ from equation 
(|2J|. This yields the results in Table- 1 which show that the 
energy cascade rate is reduced to less than 1 % of its original 
value by the time the cascade reaches k±pi ~ 32. The plasma 
parameters used to compute the damping rates in Figure- 1 are 
typical of high-speed solar wind streams at 1 AU. For these 
parameters, the electron gyro-radius occurs when k±pe ~ 1 
or, equivalently, k±pi ~ 60, and the electron inertial scale 
occurs when k±c/ujpe ~ 1 or, equivalently, k±pi ~ 30. This 
simple calculation indicates that the KAW cascade is likely 
to be dissipated before reaching electron scales. As we show 
below, a more detailed model yields results that are in good 
agreement with those in Table- 1. 

3. MODEL BASED ON CONSERVATION OF ENERGY 

The energy spectrum E{k) is defined such that E{k) dk is 
the energy per unit mass contained in the fluctuations for all 
wavenumbers between k and k + dk. The total energy con- 
tained in the entire spectrum is 



E{k)dk. 



(3) 



The cascade rate e{k) is the average energy per unit mass that 
passes through wavenumber k per unit time. The energy flows 
from low to high wavenumbers. Consider a small interval 
from kXok + dk'm wavenumber space as illustrated in Figure 
2. Energy is flowing into the boundary k from the left and out 
of the boundary k + dk from the right. Conservation of energy 
requires that the energy flowing out per unit time is equal to 





2rE(k)dk 


s(k) + — dk 

OK 


> 

Energy flowing in 
per unit time 


Energy dissipated 
per unit time 


> — 

Energy flowing out 
per unit time 



k k+dk 

Fig. 2. — For an infinitesimal interval from A: to fe + dk, the conservation of 
energy in wavenumber space requires that the energy flowing into the layer 
at wavenumber k minus the energy dissipated in the layer equals the energy 
flowing out at A; + dk. This yields equation (sj. 



the energy flowing in per unit time minus the rate at which 
energy is dissipated in the layer. Hence, 



e{k + dk) - e(k) = 2jE{k) dk, 



(4) 



where 7 < is the linear damping rate obtained from the 
KAW dispersion relation. This yields 



dk 



= 2-fE{k) 



(5) 



which expresses the conservation of energy in the cascade 
process. Note that if 7 = 0, then e = const and the energy 
cascade rate is independent of k. 

Equation Q gives one relation between e and E. A second 
is Kolmogorov's relation 



e(fc) 



kE{k) 



(6) 



where t is the energy cascade time. Komlogorov's relation 
assumes there is no dissipation of energy during the cascade 
process so that all the energy at scale k cascades to higher 
wavenumbers. When dissipation is present, a fraction of 
the total energy at scale k cascades to higher wavenumbers 
and a fraction 1 — is dissipated. The energy that is dissi- 
pated does not cascade to higher wavenumbers. In this case, 
the relation (|6]l is replaced by 



a^kE{k) 



(7) 



where (fc) is the energy attenuation factor in time t. 

It is clear from physical considerations that in the dissipa- 
tion range the energy cascade rate e must depend on the dis- 
sipation rate. Hence, Kolmogorov's relation ^ that was ini- 
tially postulated for the inertial range needs to be modified. 
What is not clear is what the precise functional form should 
be. Justification for the form (|7]l can be obtained by consider- 
ing a simple shell model. Let fc„ = 2"/co denote a sequence 
of shells in wavenumber space, where n is an integer and k^ 
is an arbitrary wavenumber The energy contained in the in- 
terval between fc„ and fc„+i is knE{kn). In the absence of 
dissipation, this energy cascades to scale kn+i in one cascade 
time T{kn) so that the energy cascade rate is 



knE(kn) 



(8) 



When there is no dissipation, e„ = const and the energy cas- 
cade rate is independent of scale. When dissipation is present. 
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the energy at scale fc„ is partly dissipated before it can cascade 
to the next level. If the energy is attenuated by a factor o? in 
time r„, then the energy cascade rate becomes 



(9) 



where a may depend on fc„. This is another way to derive 
equation (Q. 

In general, equations (|6]l and (|7|i should also contain a con- 
stant coefficient A which is dimensionless. In hydrodynamic 
turbulence it follows from dimensional analysis that this co- 
efficient is of order unity. In plasma turbulence there are other 
parameters in the problem which may cause this coefficient to 
differ from unity. For example, in the inertial range, one such 
parameter is the normalized cross-helicity Uc- Thus, equation 
(|7]l takes the final form 



e{k)^A- 



^kE{k) 



(10) 



where in many cases of practical interest the dimensionless 
constant A is of order unity. 

3.1. Complete model 

Assuming that critical balance holds for the KAW cascade, 
then the energy cascade time is equal to the wave period, r = 
2tt/lu, and the energy cascade rate (fTOl i has the form 



e{k) 



Lu{k) 
271 



Aexp{ATrj/uj)kE{k), 



(11) 



where a = cxp{2TT'y / uj) is the attenuation factor for the wave 
amplitudes in one wave period and = exp^AirjJu!) is the 
attenuation factor for the energy. When equation ( flTT l is sub- 
stituted into the conservation law Q, it follows that 
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A LO 



exp(— 47r7/aj) 



or, equivalently, 



d(loge) 1 47r7 
d{\ogk) A u 



exp(— 4717/0;). 



(12) 



(13) 



When the exponential factor on the right-hand si de is omitted, 
this equ ation is almost identical to e quation (8) in Howes et al.l 
(I2008ai) with the source term S in iHowes et al.l (12008 a') set 
equal to zero. The one minor difference is that the damp- 
ing rates 7 and u) in ( fT3] l are obtained from the hot plasma 
dispersion relation ( Stixiil992i) whereas [Howes et al.l (l2008ah 
employ the damping rates obtained from gyro-kinetic theory. 

Because w > and 7 < 0, the exponential factor in ( ITST i 
causes e{k) to decrease more rapidly than it would if the ex- 
ponential factor were omitted. When |47r7/w| > 1, the expo- 
nential drives e toward zero at a faster than exponential rate 
as can be seen in the solutions presented below. At this point, 
the dissipation becomes so effective that the energy cascade 
is abruptly terminated. However, even if the exponential fac- 
tor is omitted from equation (fTST l which is equivalent to using 
equation (|6]l in place of (|7]i, the wavenumber where the cas- 
cade rate e goes to zero is still roughly of the same order of 
magnitude. Therefore, the presence of this exponential factor 
is not crucial for the conclusions reached in this study. 

If 7(fc) and a;(fc) are known from the hot plasma disper- 
sion relation, then equation (fTsT l may be solved to find how e 



depends on k. Once the solution for £{k) is known, the en- 
ergy spectrum E{k) is obtained from equation ( fTTT i. In gen- 
eral, the ratio 7/0; depends on the angle of propagation of the 
waves which needs to be taken into account. However, when 
fc± ^ fc|| or, equivalently, when the angle of propagation is 
sufficiently close to 7r/2, then the ratio 7/0; becomes approx- 
imately independent of angle. This property is used to derive 
an approximate analytic solution below. In cases where the 
angle is not sufficiently close to 7r/2 the angle dependence 
may be taken into account as follows. Assume that the energy 
in the wavevector spectrum is concentrated near the critical 
balance curve so that the angle of propagation 6* at a given 
scale fc^ = fc is determined by the critical balance relation 



cot{e) = 



■w 



Sv I 



(14) 



,ph 



where Sv± is the electron velocity perturbation of the wave in 
the plane perpendicular to Bo and W||,ph is the parallel phase 
speed of the KAWs. Note that the parallel phase speed can be 
much larger than the Alfven speed in the kinetic regime and 
this needs to be included in the critical balance relation (fT4l i. 

According to the critical balance hypothesis, the longitudi- 
nal crossing time of two wavepackets r|| = A|| / W|| g is equal to 
the nonlinear eddy turnover time tj_ — X±/ 5v±, where W|| g is 
the parallel group speed. The parallel phase and group veloc- 
ities are approximately equal as can be seen from the approx- 
imate dispersion relation u oc k ^^ VAy/^ + {k±PiY (see, fo r 
example. iHoUwe g 1999; Cranmer & van Ballegooiienll2003h . 
In addition to equation (fT4l l. the self consistent calculation of 
the propagation angle 6 requires the relation 



kE{k) ~ ((56 J 



(15) 



where the magnetic field perturbation &b±^ is measured in ve- 
locity units. For the plasma parameters of interest here, the 
energy density of KAWs is dominated by magnetic and ther- 
mal energy (i.e., density fluctuations) with an approximate 
equipartition between magnetic and thermal energy as dis- 
cussed, for example, by Terrv et al. (2001). These two con- 
tributions have been lumped together into a single term in 
( fTSl l. Critical balance tightly couples the magnitude of the 
energy spectrum to the propagation angle through equations 
(O and ( fTSl l. For KAWs, the ratio 5v^_/Sb±^ may be deter- 
mined from the hot plasma dispersion relation using the rela- 
tions Je = —iLOEoXe ' E ~ —euQVe, whcrc Je is the elec- 
tron c urrent density and Xe is the electron susceptibility dStixl 
II992I) . This completes the specification of the model which 
consists of equations ( fTTT ). ( [13] ). ( fT4l i. and ( fTSl ). 

3.2. Analytic model 

The above model can be significantly simplified when 
k± 3> fc|| so that the ratio 7/0; becomes approximately in- 
dependent of the propagation angle 9. Simulations of in- 
compressible MHD turbulence have shown that the inequal- 
ity fcj^ ^ /c|| is usually well satisfied at the smallest inertial 
range scales and, therefore, the assumption k± ^ k^^ is justi- 
fied in the kinetic regime. It is useful to ado pt the a pproximate 
expressions for the ratio j/uj derived bv IHowes et alj (l2006l) 
using gyrokinetic t heory. Using equations (62) and (63) in 
IHowes et al.l(l2006l) one obtains 
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where a; > and a is a constant given by 
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(17) 

Here, is the electron mass, rrii is the proton mass, Tg is the 
equilibrium electron temperature, Ti is the equilibrium proton 
temperature, b = 1 + T^/Ti, and = nQkBTi/{B^/2fiQ)'^, 
where ks is Boltzmann's constant, uq is the equilibrium par- 
ticle number density, Bqis the equilibrium magnetic field, and 
fif) is the permeability of free space (SI units). Equation ( fTSI l 
is valid when k±pi ^ 1 and k±pe ^ 1. Using typical plasma 
parameters for the solar wind at 1 AU, equation (fTST i was com- 
pared to the ratio j/oj obtained from the hot plasma dispersion 
relation and found to be an excellent approximation except 
when k±pi ~ 1 where (fTSI l underestimated the damping al- 
though it remained accurate to within a factor of 2 or so. 

The substitution of equation (fTSI l into (fT2] i yields the equa- 
tion 



dk 



exp(47raA:), 



where k — k±pi. The solution is given by 



— = exp 

So 
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^4:7rakQ \ 



(18) 
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If the exponential factor is omitted from equation il3[ or, 
equivalently, equation (|6]l is used in place of (Q, then instead 
of equation (fTSl l one obtains 



d(log 



dk 



Ana 



which has the solution 



e 



exp 



Ana , 



'-(k-ko) 



(20) 
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Equations ( fT9] l and (ISTT l describe how the energy cascade rate 
changes as a result of collisionless damping. The associated 
energy spectrum can be obtained from equation ( fTTT i if the 
propagation angle d{k) is known since the propagation angle 
is necessary to compute the wave frequency uj{k). Therefore, 
the energy spectrum can only be obtained by solving the com- 
plete model consisting of equations (fTTT i. (fl3T l. ( fT4] i. and (fTsT l. 
This is not attempted here because the energy spectrum is not 
needed for the purpose of this study. 

The theoretical model developed here may also be adapted 
to study hydrodynamic turbulence. In hydrodynamics, the 
damping rate is 7 = —vk^, where v is the kinematic viscos- 
ity, and the cascade time is the nonlinear eddy turnover time 
r — 1/kv, where v = {kEy/^. The model is based on equa- 
tion (Is) and either (|6]l or (Q with a = Gxp(— i^fc^r). An ana- 
lytic solution is only possible when (|6]l is used. Calculations 
of the energy cascade rate versus wavenumber obtained using 
these hydrodynamic models show that the energy cascade rate 
appoaches zero near the Kolmogorov scale 77 ^ (ly^/eY^^. 
Moreover, the model calculations in the range ki] < 1 aie in 
reasonable agreement with the cascade rate e{k) obtained us- 
ing the model spectrum in Pope (2000). In fact, the wavenum- 
ber fc?7 ~ 1 where the energy cascade terminates in Pope's 
solution lies between the cutoff wavenumbers ki] ~ 0.8 and 
kr] ~ 1.3 for the two solutions obtained using either ^ or 
0. 
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Fig. 3. — The energy cascade rate e(fc ) of the KAW cascade in the dissipa- 
tion range obtained from equation )19t for conditions typical of high-speed 
solar wind at 1 AU: ft = 1, /3e = 1/2, and J^th.i/c = 2 X 10^*. Solutions 
are shown for A = I (blue), A = A (dashed), and A = 1/4 (dot-dashed). 
The same plot is displayed using a linear scale for the vertical axis in the 
upper panel and a logarithmic scale in the lower panel. 



4. RESULTS 

Results are now presented for typical high-speed solar wind 
in the ecliptic plane near 1 AU with the plasma parame- 
ters Vs„ > 600 km/s, Pi = 1, /?e = 1/2, and wth,p/c ~ 
2 X 10~^. This choice of pa rameters is b ased on some of our 
own work, both published (lPodestall20p9|) and unpublished , 
and results published in the literatur e dFeldman etaH [19771: 
iNewburv et al.l [19981 ISchwennll200 6'). For these parameters 
and for highly oblique angles of propagation (80° < 9 < 
90°), the damping of KAWs is negligible for small wavenum- 
bers, k±pi <C 1. KAW damping starts to become significant 
around k±pi ~ 1/2. The initial wavenumber in the model 
calculations is k±pi ~ 1/2. The results for the two differ- 
ent theoretical models, equations iT% and (|2T]) . are shown in 
Figures [3] and|4] 

The energy cascade rate e shown in Figures [3] and ID de- 
creases to zero around k±pi ~ 10 and k±pi ~ 25, re- 
spectively. When e/eq <C 1 the energy cascade terminates. 
Even though the two models yield noticeably different results, 
the wavenumber where the cascade rate becomes negligibly 
small, e/so <C 1, is of the same order of magnitude in both 
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Fig. 4. — The energy cascade rate effc ) of the KAW cascade in the dissipa- 
tion range obtained from equation )2H for conditions typical of high-speed 
solar wind at 1 AU: ft = 1, ft = 1/2, and v^^-Jc = 2 X 10"''. Solutions 
are shown for A = 1 (green), A = 4 (dashed), and A = 1/4 (dot-dashed). 
The same plot is displayed using a linear scale for the vertical axis in the 
upper panel and a logarithmic scale in the lower panel. 



Figures [3] and |4] Thus, the two models are roughly consistent 
with each other and show that the KAW energy cascade can- 
not continue beyond the wavenumber k^pi ^ 25 because the 
energy flux is almost completely damped at that point. This 
is the central conclusion of this study. The behavior of the 
solutions when the parameter A is varied are also shown in 
Figures [3] and 2] The model ( fT9] l that includes the factor 
in Kolmogorov's relation ^ is less sensitive to variations in 
the parameter A than the model ( |2TI ) that is based on Kol- 
mogorov's relation (|6]l. Taken together, the results in Figures 
[3] and |4] suggest that the KAW cascade will be completely 
damped by coUisionless Landau damping before the energy 
can reach the scale of the electron gyro-radius; the electron 
gyro-radius occurs at k^pi ~ 60 for the plasma parameters 
used here. 

The wavenumber where the KAW cascade terminates can 
be estimated from equations ( fT9l ) and ( 1211 1. The wavenumber 
where e/eg reaches some small value, say e/eq = S with 
S = 10"^, can be determined as follows. From equation (fT9] l, 
assuming Ankoa ^ 1, the wavenumber where e/eq = 6 is 
given by 



where a is defined in equation ( fTTI i. For ^ = 1 and 5 — 10"^, 
this yields k±pi ~ 10. From the second solution ( 1211 1. the 
wavenumber where e/eq = 5 is given by 



k±Pi 



(23) 



k±Pi 



— log[l + Alog(l/J)], 



(22) 



For A 1 and 6 = W^'^, this yields k±pi ~ 28. These 
formulas are very convenient for estimating the wavenumber 
where the KAW cascade terminates. 

The termination point can also be expressed in terms of the 
ratio j/u!. Using equation ( fTSI l. the condition (l22l i with A = 
1 and S = 10"^ is equivalent to \j/Ld\ = 0.14. Similarly, 
the condition ( |23] | with A ^ 1 and 5 = 10^^ is equivalent 
to \j/uj\ — 0.37. The application of theoretical models like 
those in Section 3 to hydrodynamic turbulence suggests that 
the actual termination point is somewhere between these two 
model predictions, \'y/uj\ = 0.14 and \ j/uj\ = 0.37. Hence, it 
is reasonable to conclude that the cutoff for the KAW cascade 
occurs approximately when \ j/oj\ ^ 0.25. Note that for the 
heuristic model in section 2 the cascade terminates as soon as 
= exp(47r7/w) ^ 1 which implies \'y/uj\ ^ 0.18. Thus, 
the cutoffs for the KAW cascade predicted by the heuristic 
model in section 2 and the more elaborate models presented 
in section 3 agree to within a factor of 2. 

5. CONCLUSIONS 

In this study, two different methods were employed to cal- 
culate the coUisionless damping of the KAW cascade. The 
first method demonstrates the effect of compounding on wave 
dissipation during the energy cascade process as illustrated by 
the calculation in Table- 1. For nearly perpendicular propagat- 
ing KAWs under typical solar wind conditions at 1 AU, the 
ratio 1 7/0; I is usually small from k±pi = 1 to k±pi = 10 
so that \ j/oj\ <C 1. Nevertheless, because the wave damping 
is compounded at each stage of the energy cascade process 
the energy cascade is damped more rapidly than might be ex- 
pected from an examination of the wavenumber dependence 
of 7/w. The simple heuristic model used to quantify this 
compounding effect assumes that all the energy is carried by 
KAWs, that the cascade time is equal to the wave period, and 
that damping of these waves is governed by the linear Vlasov- 
Maxwell dispersion relation. The second method based on 
equation (T3[ yields quantitatively similar results eventhough 
the precise form of Kolmogorov's relation is unknown. Thus, 
both methods support the main conclusion of this work, that 
a cascade consisting solely of KAWs cannot reach the elec- 
tron gyro-scale in the solar wind at 1 AU. This conclusion is 
also suppor ted by the somewha t different cascade model de- 
veloped by Howes et all (l2008al) . 

The conclusions are, of course, subject to some uncertainty. 
For example, if the cascade time is reduced from one wave 
period to half a wave period, then the wavenumber where the 
KAW cascade terminates is roughly twice as large. And, if 
the cascade time is reduced to one quarter of a wave period, 
then the wavenumber where the KAW cascade terminates is 
roughly four times as large. Consequently, the results are sen- 
sitive to the cascade time which is not known precisely. 

A simple expression for the wavenumber where the KAW 
cascade terminates has been obtained from the analytical so- 
lutions ( fT9l ) and ( 1211 1. For A — 1, the cutoff occurs approxi- 
mately when \ j/uj\ ~ 0.25. Hence, one must be careful when 
using damping rates obtained from the linear Vlasov-Maxwell 
wave theory to estimate the stopping point of the cascade. As 
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a consequence of the effects of compounding, the wavenum- 
ber where the cascade terminates occurs not when — 1 

but approximately when |7/a;| ~ 0.25. 

The conclusion reached in the study of ISahraoui et al.l 
(120091) should be reexamined in light of this result. 
ISahraoui et al.l (|2009) have suggested that high-frequency 
magnetic field spectra in the solar wind may be caused by a 
KA W cascade from proton to electron scales. To support this 
idea ISahraoui et alj (120091) compute dispersion curves from 
the hot plasma dispersion relation using the solar wind param- 
eters in their measurements and show that the ratio |7/w| in- 
creases to values of order \"i/oj\ — 1, indicating strong damp- 
ing, when the wavenumber reaches the scale of the thermal 
electron gyro -radius kpe = 1. Th e physical interpretation 
suggested by 'S ahraoui et al.l (|2009() appears to be untenable 
based on the model calculations presented h ere. Using the 
dispersion curves in Figure 5 of the paper by ISahraoui et al.l 
(12009) . the condition for the termination of the KAW cascade 
\"f/oj\ — 0.25 occurs when kpp ~ 30. Hence, the KAW cas- 
cade cannot reach the electron gyro-scale at kpp ~ 95 and 
cannot by itself acc ount for the h igh frequency solar wind 
spectra reported by Sahraoui et alJ (I2009D . This means that 
at least part of the high-frequency power-law sp ectrum mea- 
sured in the solar win d by Sahraoui et al. (200^. iKiyani et aT] 
( l2009i) . and lAlexandrova et al.. (.2009.) must be supported by 
some other kind of wave modes. 

The theoretical model developed here depends crucially on 
the wave damping rates computed from the linear dispersion 
relation and on the hypothesis of critical balance. Particle dis- 
tribution functions in the solar wind are usually anisotropic 
with high energy tails and sometimes other features which 
can deviate significantly from the isotropic Maxwell distri- 
butions considered here. The effects of these features on the 
wave damping rates need to be taken into account in a more 
thorough analysis. If the hypothesis of critical balance fails 
to hold or is modified somehow in the dissipation range, then 
the theory presented here would be incorrect and would have 
to be revised. A separate effect that needs to be taken into 
account for the interpretation of solar wind measurements is 
the breakdown of Taylor's hypothesis caused by the increas- 
ing phase velocities of the waves for k^pi > 1. This should 
reduce the magnitude of the spectral exponent somewhat by 
systematically shifting power to higher wavenumbers, thereby 
creating a shallower spectrum. 

It is of interest to apply the model developed here to study 
the damping of MHD turbulence in the solar corona. Be- 
cause the plasma parameters vary significantly between the 
low and high corona and the particle distribution functions 
are nonthermal, a survey of parameter space is required. This 
is an important avenue of investigation for future work. It 
is also of some practical importance since NASA plans to 
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launch a spacecraft into the high corona in the coming decade 
(iMcComas et al.1 120071) . The Solar Probe will be the first 
spacecraft to reach perihelion in the ecliptic plane at a he- 
liocentric distance of approximately 10 solar radii where the 
plasma beta is around 1/10. 

For the representative coronal parameters (3p — (3^ = 0.1 
and Tp ~ 10^ K the prediction for an isothermal electron- 
proton plasma is that the KAW cascade terminates when 
kj_pi ~ 7. For I3p = l3e = IQ-^ and Tp = 10*^ K the theory 
predicts that the KAW cascade terminates when k±pi ^ 2.5. 
It should be noted that for these plasma parameters the ap- 
proximation ( fTSI l overestimates the dampins rate at k±pi ^ 1 




Fig. 5. — The ratio ■y/uj obtained from the hot plasma dispersion relation 
for a propagation angle 9 = 89.9 degrees is compared to the approximation 
(16) used in the theoretical model (dashed) for fip = Pe = 10~^ and Tp = 
K (red), /3p = = 10~^ and Tp = lO"^ K (blue), and /3p = 1, 
/3e = 1/2 and «^th,p = 60 km/s (green) 



by a factor of ^ 1.5 and, consequently, these preliminary esti- 
mates could be too small by a factor of 2 or so. Nevertheless, 
these preliminary calculations suggest that the KAW cascade 
cannot reach the electron gyro-scale in the solar corona. 

The cascade terminates more rapidly at low j3p because for 
a fixed value of fcj^Pi the ratio |7/a;| increases as Pp decreases. 
This is shown in Figure-|5] where the ratio 7/0; obtained from 
the hot plasma dispersion relation is compared to the approx- 
imation ( fT6b used in the theoretical model. Note from equa- 
tions ( fT6] l and ( fTTj i that when — Tp and (3p <s: 1, if k±p.^ is 

held fixed {k±pi ^ 1), then \ j/uj\ (x /3p even though for 
a fixed value of k±c/ujpi the ratio |7/a;| is independent of /3p. 

This work was supported by the NASA Solar and He- 
liospheric Physics Program, the NASA Heliospheric Guest- 
Investigator Program, the NSF SHINE Program, and by the 
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